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We report lattice-volume independence of low moments of nucleon structure functions from the 
coarse RIKEN-BNL-Columbia (RBC) and UKQCD joint dynamical (2+l)-flavor domain-wall 
fermions (DWF) ensembles at the lattice cut off of a~ 1 ~ 1.7 GeV. The isovector quark momen- 
tum fraction, (x) u -d> ar, d helicity fraction, (x)ah-a</> both fully non-perturbatively renormalized 
are studied on two spatial volumes of (~ 2.7fm) 3 and (~ 1.8fm) 3 . Their naturally renormalized 
ratio, (x) u -d/ (x}au-M< is not affected by any finite-size effect. It does not depend strongly on 
light quark mass and does agree well with the experiment. The respective absolute values, fully 
non-perturbatively renormalized, do not show any finite-size effect either. They show trending 
toward the respective experimental values at the lightest up- and down-quark mass. This trending 
down to the experimental values appears to be a real physical effect driven by lighter quarks. The 
observations are in contrast to the huge finite-size effect seen in the axial-current form factors. 
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1. Introduction 

We report lattice-volume independence of low moments of nucleon structure functions. In 
Lattice 2008 we reported our lattice numerical calculations of some low moments of nucleon struc- 
ture functions [1] with the RBC and UKQCD joint dynamical (2+l)-flavor domain-wall fermions 
(DWF) ensembles at the lattice cut off of cT l ~ 1.7 GeV [2]. In contrast to companion form factor 
calculations [3, 4], we did not see any sign of the huge finite-size effect in the structure functions. 
However those calculations then were complete with only the larger, (2.7fm) 3 volume. Now we 
have completed the calculation with the smaller, (1.8fm) 3 volume, and are more confident that the 
structure functions do not suffer from the finite-size effect so much as the form factors. 

The nucleon structure functions are measured in deep inelastic scattering of leptons off nu- 
cleons [5, 6], the cross section of which is factorized in terms of leptonic and hadronic tensors, 
1^ V W^ V . Since the leptonic tensor, l^ v , is known, the cross section provides us with structure 
information about the target nucleon through the hadronic tensor, W^ v , which is decomposed into 
symmetric unpolarized and antisymmetric polarized parts: 

wW(x,Q 2 ) = ieW°q p (^Ui(*,<2 2 ) + g 2 (x,Q 2 ))-^g 2 (x,Q 2 )^ , (1.2) 

with kinematic variables defined as v = q-P, S 2 = —M 2 , and* = Q 2 /2v, and Q 2 = \q 2 \. The unpo- 
larized structure functions are F\(x, Q 2 ) and F 2 (x,Q 2 ), and the polarized, gi(x,Q 2 ) and g 2 (x, Q 2 ). 
Their moments are described in terms of Wilson's operator product expansion: 

2 f dxx^F^Q 2 ) = £ cS(y%( J u) + 0(l/G 2 ), 

J0 q=u,d 

f dxx"- 2 F 2 {x,Q 2 ) = £ cK(*'%Gu) + 0(l/<2 2 ), 
2 f dx^ gl (x,Q 2 ) = £ e { ?l(x») Aq (ii) + 0(\/Q 2 ), 

+ <9(l/<2 2 ), (1.3) 

' q=u,a 

where we suppressed the (jU 2 /Q 2 ,g(jU)) dependence of the perturbatively known Wilson coeffi- 
cients, a, c 2 , e\, and e 2 . The moments, {x") q (lJ,), (Oa<?(a0 an< ^ ^»(a0 we calculable on the lattice 
as forward nucleon matrix elements of certain local operators. 

In the following we report our lattice numerical calculations of two lowest isovector mo- 
ments of the structure functions, the quark momentum fraction, {x) u -d{pL), the helicity fraction, 
(x)Au~Ad(l^), and their ratio. We note the results for the tensor charge, (1) Su-Sdil 1 ), which probes 
transverse spin structure of nucleon [7, 8], and twist-3 coefficient d\ of the g 2 polarized structure 
function, were reported in Lattice 2008 [1]. 
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2. Formulation 

We use the standard proton operator, B = e a bc{u^Cy^db)u c to create and annihilate proton 
states. We Gaussian-smear this operator for better overlap with the ground state with both zero 
and finite momentum, as we conducted a companion form-factor calculations [3, 4]. Since the up 
and down quark mass are degenerate in these calculations, the isospin symmetry is exact. This is 
of course a well-known good approximation. We project the positive -parity ground state, so our 
two-point proton function takes the form 

C2pt(0 = [ — (Bp{hmk)B a (/source)), (2-1) 
a,p V 1 J a/3 

with t = ^sink — Source- We insert an appropriate observable operator 0(q, t') at time t' = ? SO urce + T, 
< T < t, and possibly finite momentum transfer q, to obtain a form factor or structure function 
moment three-point function, 

c \p.{t^A) = J^r a p(Bp(t s i nk )0(q,t')B a (t somce )}, (2.2) 

1+7; 1+7? 
with appropriate projection, F = — - — , for the spin-unpolarized, and F = — - — iYsYk,k / 4, for 

the polarized. Ratios of these two- and three-point functions give plateaux for < % < t that 

C 3 r f (f,T) 

give the bare lattice matrix elements of desired observables: e.g. (0) bare = 3pt at q 2 = 0. 

C2pt(0 

We renormalize the structure function moments by Rome-Southampton RI-MOM non-perturbative 
renomalization prescription [9]. The good continuum- like flavor and chiral symmetries of domain- 
wall fermions [10, 1 1, 12] are very useful here in eliminating unwanted lattice-artifact mixings that 
are present in many other fermion schemes. See our earlier publications, ref. [13], and references 
cited there in for further details. 

The RBC-UKQCD joint (2+l)-flavor dynamical DWF coarse ensembles [2] are used for the 
calculations. These ensembles are generated with Iwasaki action at the coupling j3 =2.13 which 
corresponds to the lattice cut off of about a -1 = 1.73(3) GeV. There are two lattice volumes, 
16 3 x 32 and 24 3 x 64, corresponding to linear spatial extent of about 1.8 and 2.7 fm. The dynamical 
strange and up and down quarks are described by DWF actions with the fifth-dimensional mass of 
M5 = 1.8. The strange mass is set at 0.04 in lattice unit and turned out to be about twelve percent 
heavier than physical including the additive correction of the residual mass, m res = 0.003. The 
degenerate up and down mass is varied at 0.03, 0.02, 0.01 and 0.005. We summarize the gauge 
configurations and pion and nucleon mass of each ensemble in Table 1 . 

There are two important sources of systematic error: finite spatial size of the lattice and excited 
states contamination: Chiral-perturbation-inspired analysis of the former for meson observables 
suggests the dimensionless product, m n L, of the calculated pion mass, m % and lattice linear spatial 
extent, L, should be set greater than 4 to drive the finite-volume collection below one percent, 
and the available lattice calculations seem to support this. While our present parameters satisfy this 
condition, it should be emphasized that this criterion is not known sufficient for baryon observables. 
Indeed our companion calculations on nucleon isovector form factors [3, 4] strongly suggests much 
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nifa 


# of config.'s 


meas. interval 


N sources 


m n (GeV) 


niN (GeV) 


0.005 


646+286 


10 


4 


0.3294(13) 


1.154(7) 


0.01 


356 


10 


4 


0.4164(12) 


1.216(7) 


0.02 


98 


20 


4 


0.5550(12) 


1.381(12) 


0.03 


106 


20 


4 


0.6681(15) 


1.546(12) 



Table 1: Number of gauge configurations and pion and nucleon mass. 



larger volumes with m n L ~ 6 — 7 are necessary. Here we test this by performing structure function 
calculations on two lattice volumes, (~ 2.7fm) 3 and (~ 1.8fm) 3 

One should adjust the time separation, t, between the nucleon source and sink appropriately so 
the resultant nucleon observables are free of contamination from excited states. The separation has 
to be made longer as we set the quark masses lighter. Here our previous study with two dynamical 
flavors of DWF quarks [13] help as the lattice cutoff is similar at 1.7 GeV: By comparing the 
isovector quark momentum fraction, at the up/down mass of 0.02, where pion mass is 

about 500 MeV and nucleon 1.3 GeV, a clear systematic difference is seen between the shorter time 
separation of 10, or about 1.16 fm, and longer 12, or 1.39 fm: the former averages to 0.236(9) while 
the latter 0.195(17). Although the former appears statistically better, as the two-point function that 
provides the normalization has not decayed much yet, it is not manifestly free of excited-state 
contamination that is expected more at shorter time separation. The latter suffers statistically, 
as the two-point function that provides the normalization has decayed, but is certainly freer of 
excited-state contamination than the former: Of the two, we have to choose the latter. If larger 
statistics is required, then we have to collect as much statistics as necessary for the latter. Ideally 
we should collect more statistics at even larger source-sink separation to confirm the excited-state 
contamination has been eliminated by the separation of 12. 

Since in the present ensembles we explore much lighter up/down quark mass than this two- 
flavor example, with pion mass as light as 330 MeV and nucleon 1.15 GeV, the source/sink time 
separation must be chosen longer: longer than 1 .2 fm at least. For the lightest up/down quark mass, 
rrifa = 0.005, in the present (2+l)-flavor dynamical DWF calculation, the nucleon signal begins to 
decay at t = 12, or about 1.4 fm [14] (see Figure 1): hence it is about the longest distance we can 
choose without losing the signal, and hence about as free of excited-state contamination as we can 
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Figure 1: Nucleon effective mass at quark mass nifa = 0.005 in the present calculation with (2+l)-flavor 
dynamical DWF quarks. The nucleon signal begins to decay at time separation 12. By which time we are as 
free of excited state as possible while maintaining the nucleon ground state. 
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Figure 2: Ratio, (x) u - d /{x) Au-Mj °f isovector momentum and helicity fractions, naturally renormalized 
on the lattice, from (2.7fm) 3 (x) (1.8fm) 3 (+), plotted against pion mass squared in GeV 2 . They are in 
excellent agreement with the experiment which is denoted by a square (□). No discernible dependence on 
volume nor pion mass can be detected. 

go. We found the bare three-point function signals for quark momentum and helicity fractions for 
this source-sink separation of t = 12 are acceptable [1]. Thus we choose the source-sink separation 
of 12, or about 1.4 fm for the present work, as the two-flavor example suggests strongly that any- 
thing shorter would suffer from excited-state contamination. To confirm this choice is sufficiently 
long to eliminate the excited-state contamination remains a future problem. 

3. Results 

The ratio, (jc) m _<//(jc)am-a<*> of the isovector quark momentum fraction to the helicity fraction, 
is naturally renormalized in DWF lattice QCD because the two quantities share a common renor- 
malization [15]: the momentum fraction, (x) u -a, the first moment of the unpolarized structure 
functions, and helicity fraction, (jc)A»-Arf> the first moment of the g\ polarized structure function, 
are related by a chiral rotation and the DWF preserve the chiral symmetry. Thus the ratio calcu- 
lated on the lattice is directly comparable with the experiment. Our results for this ratio do not 
show any discernible dependence on the up/down quark mass, and are in excellent agreement with 
the experiment (see Figure 2.) This is in contrast to the similarly naturally renormalized ratio of 
gA/gv °f elastic form factors which at the lightest point deviates significantly from heavier mass 
results and the experiment [3]. This suggests individual moments of inelastic structure functions 
such as the momentum fraction, (x) u -d, and helicity fraction, (x)a u -m, may not suffer so severely 
from the finite-size effect that plagues elastic form factor calculations. 

The absolute value of the isovector quark momentum fraction, (*)„_</, fully non-perturbatively 
renormalized and run to MS at 2 GeV, Z MS (2GeV) = 1.15(4), is shown in Figure 3. On the larger, 
2.7 fm, volume, the three heavier points stay roughly the constant which is about 70 % higher at 
~ 0.26 than the experiment, about 0.15. This behavior is not so different from old RBC quenched 
results [15] wtih similar up/down quark mass. The lightest point shows a sign of deviation away 
from this heavy constant behavior, but in contrast to the form factor deviations, is trending toward 
the experiment. Since lighter quark can more easily share its momentum with other degrees of 
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Figure 3: Absolute value of isovector quark momentum fraction, (jc) u _</, fully non-perturbatively renor- 
malized. Plotted against the pion mass squared in GeV 2 (left) in comparison with the experiment which is 
denoted by a square (□), and against m K L (right). Results from the larger, L = 2.7 fm volume are described 
by x symbols and the smaller L = 1 .8 fm volume by +. 



freedom, this trending toward the experiment may well be a real physical effect: It does not have to 
be affected by the finite spatial size of the lattice. This question is now clarified by the calculations 
on the small, (1.8fm) 3 , volume: at rrifa = 0.01 on the smaller volume of 1.8 fm, with m n L smaller 
than that at the mja = 0.005 on 2.7 fm volume, the calculated momentum fraction is in agreement 
with the large volume result at the same mass. This absence of the finite-size effect scaling in m K L 
can also be seen in the right pane of the figure where the same quantities are plotted against this 
scaling variable. 

The isovector quark helicity fraction, (jc)ah-/W> is also fully non-perturbatively renormalized 
and run to MS at 2 GeV (See Figure 4.) The renormalization, Z MS (2GeV) = 1.15(3), agrees very 
well with the momentum fraction one. As can be expected from the constant behavior of the ratio, 



m 



0.1 0.2 0.3 0.4 



Figure 4: Absolute value of isovector quark helicity fraction, (x}au-M< fully non-perturbatively renormal- 
ized. Plotted against the pion mass squared in GeV 2 in comparison with the experiment which is denoted 
by a square (□). Results from the larger, L = 2.7 fm volume are described by x symbols and the smaller 
L= 1 .8 fm volume by +. 



6 



2+lfDWF nucleon structure function 



Shigemi Ohta 



(*)»-<*/(*) Ah-m> this denominator observable exhibits very similar behavior to the numerator one 
discussed in the above: the three heavier points stay roughly the constant which is about 70 % 
higher than the experiment. The lightest point then shows a sign of deviation away from this 
heavy constant behavior and trends toward the experiment. This trending toward the experiment is 
independent of the lattice volume and of the form factor scaling variable, m n L: It is likely a real 
physical effect arising from light quark mass. 

4. Conclusions 

A naturally renormalized observable with DWF quarks because of well-preserved chiral sym- 
metry, the ratio, (x) u -d / (x}au-AcI, of the isovector quark momentum and helicity fractions, does not 
depend strongly on the degenerate up- and down-quark mass and does agree well with the experi- 
ment. This is in contrast to a similarly naturally renormalized ratio, the axial charge, gA/gv> where 
a huge finite-size effect takes the calculation away from the experiment at ligher ud-quark mass 
values [3] . The absolute values of the fractions, {x) u -d and {x)\ u -m, fully non-perturbatively 
renormalized, show trending toward the respective experimental values at the lightest ud-quark 
mass value, away from constancy at heavier values. Comparison of the results on the two different 
lattice volumes does not reveal any finite-size effect. In particular, we rule out huge finite-size 
effect scaling in m n L seen in the axial charge [3]. Thus the trending down to the respective experi- 
mental values of these two isovector fractions appears to be a real physical effect driven by lighter 
ud-quarks, in contrast to the finite-size effect in the axial-current form factors. 

I thank J. Zanotti, T. Yamazaki, and other RBC and UKQCD members. RIKEN, BNL, US 
DOE, Edinburgh University and UK PPARC provided facilities essential for this work. 

References 

[1] S. Ohta and T. Yamazaki, [RBC+UKQCD], PoS (LATTICE 2008) 168, arXiv:08 10.0045 [hep-lat]. 

[2] C. Allton et al, [RBC+UKQCD], Phys. Rev. D 78, 1 14509 (2008), arXiv:0804.0473 [hep-lat]. 

[3] T. Yamazaki et al, [RBC+UKQCD], Phys. Rev. Lett. 100, 171602 (2008), arXiv:0801.4016 [hep-lat]. 

[4] T. Yamazaki et al, [RBC+UKQCD], Phys. Rev. D 79, 1 14505 (2009), arXiv:0904.2039 [hep-lat]. 

[5] E. D. Bloom et al, Phys. Rev. Lett. 23, 930 (1969). 

[6] M. Breidenbach et al, Phys. Rev. Lett. 23, 935 (1969). 

[7] M. Anselmino et al, Phys. Rev. D 75, 054032 (2007), arXiv:hep-ph/0701006. 

[8] RHIC Spin Collaboration, http://spin.riken.bnl.gov/rsc/report/spinplan_2008/spinplan08.pdf (2008). 

[9] C. Dawson et al, Nucl. Phys. B 514, 313 (1998), arXiv:hep-lat/9707009. 

[10] D. B. Kaplan, Phys. Lett. B 288, 342 (1992), arXiv:hep-lat/9206013. 

[11] Y. Shamir, Nucl. Phys. B 406, 90 (1993), arXiv:hep-lat/9303005. 

[12] V. Furman and Y. Shamir, Nucl. Phys. B 439, 54 (1995), arXiv:hep-lat/9405004. 

[13] H. W. Lin et al, [RBC], Phys. Rev. D 78, 014505 (2008), arXiv:0802.0863 [hep-lat]. 

[14] T. Blum, [RBC+UKQCD], PoS (LATTICE 2008) 096. 

[15] K. Orginos, T. Blum and S. Ohta, [RBC], Phys. Rev. D 73, 094503 (2006), arXiv:hep-lat/0505024. 



7 



